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The role of quantum coherence loss in mixed quantum-classical dynami- 
cal systems is explored in the context of the theory of quantum decoherence 
introduced recently by Bittner and Rossky. (J. Chem. Phys. 103, 8130 
(1995)). This theory, which is based upon the consistent histories interpreta- 
tion of quantum mechanics, introduces decoherence in the quantum subsys- 
tem by carefully considering the relevant time and length scales over which 
one must consider the effects of phase interference between alternative his- 
tories of the classical subsystem. Such alternative histories are an integral 
part of any quantum-classical computational scheme which employ transi- 
tions between discrete quantum states; consequently, the coherences between 
alternative histories have a profound effect on the transition probability be- 
tween quantum states. In this paper, we review the Bittner-Rossky theory 
and detail a computational algorithm suitable for large-scale quantum molec- 
ular dynamics simulations which implements this theory. Application of the 
algorithm towards the relaxation of a photoexcited aqueous electron compare 
well to previous estimates of the excited state survival time as well as to the 
experimental measurements. 



I. INTRODUCTION 



A key issue which has emerged in describing electronic transitions in the condensed phase 
is the proper treatment of the short-hved phase coherence of the electronic and nuclear (bath) 
wavefunction. The essence of the effect becomes clear if one considers the behavior of an 
electronic quantum system in a classical bath. In such a mixed quantum-classical system, 
the classical dynamics of the bath nuclei follow along a given adiabatic potential surface 
associated with one of the eigenstates of the electronic Hamiltonian. Thus the potential 
felt by the nuclear degrees of freedom will depend strongly upon the quantum state of the 
electronic degrees of freedom. The fundamental distinction between quantum mechanics and 
classical mechanics is that a quantum system can evolve into a coherent linear superposition 
of states. This leads to a rather profound ambiguity in describing the forces between the 
quantum and classical subsystems. Furthermore, quantum mechanics forces us to consider 
the effect of all alternative histories or pathways, including those of the classical bath, 
when computing quantum transition probabilities. Thus, as the electronic wavefunction 
evolves from an initially pure eigenstate to a coherent superposition of eigenstates, various 
alternative pathways for the bath particles begin to emerge, each associated with dynamics 
starting from an initial nuclear configuration and electronic state and ending in a different 
final nuclear configuration and final electronic eigenstate. As the nuclear dynamics for 
the different quantum (electronic) states diverge, phase coherence between the alternative 
histories is lost due to the increasing differences between the classical actions associated with 
each alternative nuclear path. This effect is known as quantum decoherence.0 

Because quantum decoherence acts to suppress the formation of superposition states, 
it diminishes profoundly the transition probability between quantum states coupled by the 
nuclear dynamics. When the loss of quantum phase coherence between the electronic states 
is neglected and no longer enters directly into the computation of making an electronic 
transition, which is a typical approximation made in mixed quantum-classical treatments, 
one can show that the result will be an incorrect estimation of the transition probabilities 
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and other associated physical observables.i Because of the tremendous utihty of mixed- 
quantum classical treatments in simulating condensed phase phenomena, it is of tantamount 
importance to be able to properly incorporate the effects of quantum decoherence into these 
simulations. 

Recently, we have presented a theory of quantum decoherence which is suitable for con- 
densed phase quantum molecular dynamics simulations.! This approach is based upon a 
interpretation of quantum mechanics introduced by Griffithsi, Omne^li, and Gell-Mann 
and Hartlei^ over roughly the past 10 years. This interpretation of quantum mechanics, 
termed "consistent" or "decoherent" histories, has generated a great deal of attention in 
the field of quantum cosmolog We have also applied and extended the formalism to 
study coherence effects in physically realistic condensed phase problems This work has 
lead to a number of advances in designing computational algorithms which incorporate a 
consistent description of quantum decoherence and to a deeper understanding of the role 
that quantum decoherence plays in condensed phase chemical processes and molecular level 
descriptions of the decoherence pro cess. These theoretical advances have also offered 
a tangible solution to a specific long standing puzzle regarding the solvent isotope effect 
on the non-radiative lifetime of an excess electron in H2O and D2O and offer a plausible 
explanation of the lack of an observed solvent isotope effect .lil^ 

In this paper, we give the details of our computational algorithm based upon consistent 
quantum histories and its application in a full scale quantum molecular dynamics simula- 
tion of the relaxation of an excited excess electron in H2O. The simulations presented in 
this paper represent the first application of the Bittner-Rossky decoherence theory in a full 
scale simulation of a condensed phase system. The aqueous electron provides perhaps the 
simplest example of electronic dynamics in the condensed phase. At the same time, aqueous 
electrons are of particular importance since they play a prominent role in the broad area 
of radiation chemistry of water. In spite of intensive experimental and theoretical study 
since identification of the hydrated electron over 30 years ago, many dynamical features of 
electron solvation have remained incompletely understood, largely due to the extremely fast 



time scales on which electronic relaxation occurs.El Newly developed femtosecond spectro- 
scopic techniques are now providing glimpses of the details of condensed phase electronic 
dynamics. From a theoretical standpoint, the aqueous electron is an ideal test case for new 
condensed phase theories. Recent theoretical advances in treating condensed phase elec- 
tronic dynamicsEi^ has led to an exceptional interplay between theory and experiment, 
helping to unravel many subtle dynamical aspects of this system. 

The remainder of this paper proceeds as follows: In the next section we review the 
consistent histories methodology as it applies toward developing a computational algorithm 
for condensed phase simulations. We then give the results of quantum MD simulations of 
the non-adiabatic relaxation dynamics of an excited excess aqueous electron which employ 
the consistent histories algorithm and compare the present results to previous estimates 
of the non-adiabatic electronic transition rate by Schwartz, Bittner, Prezhdo, and Rossky 
(SBPR)lll, to earlier results presented by Schwartz and Rossky which employed a constant 
1 fs coherence time scale throughout as well as to the experimental data of Kimura et 
a/..0Hli Finally, we comment upon the utility of the present algorithm and discuss its limi- 
tations as well as future extensions and improvements. 



II. THEORETICAL METHODS 

A. Coherence Between Switching Paths 

Let us consider the evolution of an initial quantum state \i{Ro)), taken as an adiabatic 
eigenstate of the quantum mechanical Hamiltonian for initial bath (or classical) configura- 
tion Ro, along a switching path, which begins at Ro. Along the course of this path, we 
determine at various time intervals whether or not a quantum transition has occurred ac- 
cording to a stochastic switching criteria and modifying the classical dynamics accordingly. 
The fundimental notion of having classical particles switch between discrete quantum states 
was originally pioneered by TullyS in the early 1970's. Subsequently, various computa- 
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tional schemes have been developed which incorporate quantum switching into the classical 
dynamics and improve upon Tully's original idea. Most notably are the stationary phase 
surface hopping algorithms developed by Webster, Friesner, and Rossky^S and further by 
Coker and co-workerJ^lllliiil^ as well as the Molecular Dynamics with Quantum Transitions 
algorithm of Tully and co-workersilS. The classical switching path contains a record of the 
outcome of each switching attempt and can thus be written as a time ordered sequence of 
events, 

R^{t) = {R:°,---,RT---,R:-}. (1) 

The superscript aj denotes the switching outcome at time-step j and ao = i corresponding 
to our choice of the initial quantum state. Two such paths are shown schematically in Fig. |l| 
where we plot the eigenenergy of the occupied state along the path as a function of time. 
Changes in the quantum state imply that there is a corresponding sudden change in the 
forces exerted on the classical particles over the course of its evolution. The result is that 
different sequences of switching events will lead to rapidly diverging paths. 

Along a given switching path, -R"(t), the partial transition amplitude between an initial 
and final quantum state, \i{Ro)) and \ f{Rf)) respectively, is given by 

T,f[R'^{t)] = (/(i?^)|e-'^'l'"^^'2f^"(^^^K(/2o))e+'^[^"(*)] 

= f/,/[i?'^(t)]e+'^[^°W], (2) 

where Uif[R°'{t)] is the transition amplitude for the quantum subsystem (with Hamiltonian 
Hq) and ^[/^"(t)] is the classical action computed along the switching path R"{t), i.e. 

S[R{t)] = £ dt^'^Ritr - Vcm))] . (3) 

Here, Vc is the interaction potential coupling the classical particles. The probability of 
starting in some initial state |«(to)) ending up in a given final quantum state, \f(tf)), at 
some time t/ is computed by summing over all possible switching paths which connect the 
initial and final quantum states. Averaging over initial configurations, the reduced transition 
probability for a mixed quantum-classical system is 




(tf) = ( E T,fmt)]T:,[R^{t)]p.{Ro) 



\{R'-{t),R0it)} 



(4) 



where Pi{Ro) is the probabihty density of being in the initial state, the average, (■ ■ ■)o, is 
taken over initial configurations, and the sum is over pairs of switching paths 



which start at the initial configuration, Ro, with the quantum state in the initial state and 
end at any final configuration with the quantum state in state |/) at tf. 

According to this last equation, the transition probability for the quantum subsystem is 
dependent upon interferences between the various alternative pairs of switching pathways. 
Thus, in order to correctly compute the populations of the electronic states at some later 
time, one must perform a sum over all alternative paths that the bath can take. In a 
practical sense, this amounts to launching a swarm of trajectories from each initial starting 
point and computing the switching probabilities (i.e. branching ratios for the swarm) using 
Eq. ^. However, this approach has the distinct disadvantage that one must initialize a large 
enough swarm of trajectories in order to effectively sample the distribution of switching 
paths originating from each initial configuration. So, while formally correct and tractable in 
small dimensional systems, the tremendous computational overhead required to implement 
Eq. ^ directly for a large number of classical variables is prohibitive. 

Next, we present an alternative derivation of the quantum survival and transition prob- 
abilities which uses only a single switching path explicitly. The formalism is motivated by 
the observation that at short times, there will be a significant contribution from the phase 
interferences between two alternative paths with similar histories, necessitating a coherent 
sum over paths, while at longer times the action difference between two paths will become 
very large and phase interference effects from alternative paths will become negligible. The 
time scale separating these two regimes is the quantum coherence time scale. 



{i?"(t),i?^(t)|i?"(0) = i?o,i?^(0) = i?o} 



(5) 
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B. Consistent Histories 

Let us define the space of paths, as consisting of all continuous paths connecting some 
initial system/bath state to the final state i?j labeled by a set of quantum numbers, {cxj} 
and a classical bath variable, R. Define also a set of subspaces of Z labeled by A(i?") which 
are the spaces of paths taken by the quantum subsystem which are parameterized by the 
paths taken by the bath. The full space of paths is the continuous union of all the various 
subspaces. 

Z = |+)A(i?") (6) 

a 

These paths are completely fine grained histories since the values of the each paths can be 
specified at all times and one can compute the partial amplitude along any one of the paths 
in Z. Furthermore, by integrating over the individual "quantum" subspaces one obtains 
Eq. I, 

JA{_R«) 

= [/,y[i?"(t)]e+^^[^'^Wl. (7) 

A fundamental property of quantum mechanics is that fine grained histories can not 
be assigned probabilities, only amplitudes. In order to make predictions based upon the 
theory of probabilities, coarse graining of histories must be introduced. Coarse graining is 
accomplished when a "reduced" history is constructed by summing over sets of fine grained 
histories at a particular time. Coarse grained histories can be assigned probabilities, that 
is, integrating over coarse grained histories is equivalent to integrating over probabilities. 

Coarse graining of histories arises naturally whenever we try to mix quantum with clas- 
sical dynamics. As different classical switching paths emerge from a common origin, the 
space of quantum paths parameterized by one switching path will diverge from the space 
of quantum paths parameterized by another. As the "overlap" between subspaces decays 
(corresponding to different classical dynamics), the originally fine grained space of paths, Z, 
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becomes coarse grained into separate subspaces, {A(i?°)} as time goes on. Consequently, 
transition probabilities are computed by first integrating coherently over the fine grained 
sets of paths to get the transition amplitude say along one classical switching path, then 
summing incoherently (i.e. adding as ordinary probabilities) over coarse grained sets. 

Quantum decoherence effects can be consistently incorporated into mixed quantum- 
classical systems by recognizing that restricting the quantum evolution to given classical 
pathways is equivalent to making a series of quantum measurements on the total system. 
Furthermore, the quantum coherence time is the time scale which characterizes how often 
such "measurements" occur!. The classical path sequence shown in Fig. |I| is an example of 
a quantum mechanical history. Along this history, the time evolution operator for the quan- 
tum system can be written equivalently as a time ordered sequence of alternating quantum 
projection operators and unitary evolution operators 

C[R'{t)] = UnP""-' ■ ■ ■ UiP'''Uo (8) 

where the projections at each time interval are members of complete sets representing the 
total set of possible outcomes, 

E^"' = l' (9) 

and 

Un = e -^'"-i (10) 

is the unitary evolution operator for the quantum wavefunction along a segment of the 
switching path. 

The decoherence time scale sets a characteristic time interval between subsequent appli- 
cations of the projection operators. This time scale is roughly the time scale for alternative 
stationary trajectories of the nuclear bath to diverge sufficiently such that the difference in 
action between the two paths in Eq. ^ is large. In Ref.l2 we suggested that coherence time 
scale is related to the time scale for the decay of the overlap integral between nuclear wave- 
functions evolving on different adiabatic potential energy surfaces. During the coherence 
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interval, transition amplitudes are added coherently according to the rules of quantum me- 
chanics. Application of the projection operators destroys the coherences and one is left with 
transition probabilities which are then added according to the rules of standard probability 
theory. Within this approach, Eq. ^ is equivalent to 

P^^At) = ^ E Un\C[R^{t)MRo)C^[R'mfn)) 

= ( E (/n|f>nP"'^-^---P"^f>i/5°f>|P"^---P""-^?7l|/^ , (11) 

\{B^{t)} 

where 1/^) = |/(-R"")) is the final quantum state at the end of path R^it). 

When the time intervals 5tn = tn — tn-i are members of a Poisson distribution, the 
probability of maintaining coherence over a given interval, t, follows from the exponential 
deviate 

7r(t) = / ds- = 1 - e-*/", (12) 

JO r 

and T is the characteristic decoherence timescale. Thus, over a short time interval, 5t, where 
the probability of collapsing the wavefunction is 

T,{t) = 5t/T + 0{5t'), (13) 

the reduced density matrix of the quantum system evolves as 

p{t + 5t) = (1 - 6t/T){p{t) + i5t[H,p\) + 5^/r7^[p]. (14) 

where TV^p] is a "reduction mapping" of the density matrix p which projects out the diagonal 
elements of the quantum density matrix (the populations) killing off the coherences between 
the quantum states. 



Taking the limit oi 5t ^ Q in Eq. O and identifying r as the decoherence timscale td 



we arrive at the master equation for the quantum density matrix. 



m = i[HiR),p{R)] - -l(p - nipiR)]), 

= C{R)p--{p{R)-n[p{R)]). (15) 



The first term in this equation contains the Liouvillian of the quantum system, C{R) = 
i[H{R), ]. Evolution of the quantum system under this term alone is unitary and non- 
dissipative. The second term introduces quantum decoherence into the dynamics of the 
quantum subsystem. The coherences originally in the quantum subsystem decay due to the 
series of measurements imposed by the environment. Although energy exchange does not 
appear explicitly in Eq. |15|, energy exchange between the system and the bath is explicitly 
included through dynamics of the bath variable R(t) treated separately. 

In other theories of quantum relaxation, such as the spin-boson model US, the Redfield 
equations or Liouville space methods^, both decoherence and dissipation are treated 

implicitly through effective equations of motion for the bath, thus losing the molecular level 
information about the underlying dynamics of the bath. In our treatment, dissipation is in- 
cluded explicitly through the classical molecular dynamics of the condensed phase medium, 
thus we are able to retain a molecular level description of the bath while quantum deco- 
herence is treated implicitly in order to avoid summing over alternative pairs of classical 
paths. 

It is important to note that the form of the projections and the coherence time scale 
(and the resulting master equation) are directly related to the forces coupling the quantum 
system to the bath and hence are very dependent upon the choice of basis used to represent 
the quantum subsystem. In the present application of our theory, we assume that the 
projections periodically resolve the quantum subsystem into the adiabatic states with a 
time scale determined by the average decay time of the Franck-Condon overlap between 
nuclear wavefunctions initialized on the different adiabatic surfaces at the instantaneous 
nuclear coordinates.^ A more detailed analysis of this assumption is forthcoming.0 

Perhaps the most important physical consequence of decoherence for chemical physical 
applications is the overall reduction of the transition probabilities between different quantum 
states. In order for a quantum system to make a transition from one state to another, the 
system must form a coherent superposition between the states. The ease with which these 
coherences form reflects the strength of the coupling between the two states. Decoherence, 
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on the other hand, reflects the strength of the system-bath couphng. When a quantum 
system evolves into a superposition state, the destructive interferences between divergent 
alternative switching paths leading to the final state destroys the coherence in the quantum 
subsystem, hence, reducing the likelihood of making a transition into the final state.i 



III. SIMULATIONS 

The simulation techniques employed here are nearly identical to those used in earlier 
work by Murphrey and Rossky studying both the relaxation of electrons photo-injected into 
neat water0§i as well as the present case of photoexcitation of equilibrium hydrated elec- 
trons by Schwartz and Rossky. The exceptions regard the treatment of the quantum coher- 
ences detailed below. Briefly, the model consists of 200 classical SPC water molecules 
with the addition of internal flexibility^ and one quantum electron in a cubic cell of side 
18.17 A (corresponding to a solvent density of 0.997 g/ml) with standard periodic bound- 
ary conditions at room temperature. The electron-water interactions were described with a 
pseudo-potential ,0 and the equations of motion integrated using the Verlet algorithm with 
a dynamical time step of 1 fs in the microcanonical ensemble. @ The adiabatic eigenstates at 
each time step were computed via an efficient iterative and block Lanczos scheme utilizing 
a 16^ plane wave basisi^ The lowest 6 eigenstates were used as a basis for propagating the 
quantum wavefunction and to compute the stationary phase paths for the nuclear dynamics. 
Explicit details of these methods are discussed elsewhere.@ 



A. Consistent Histories Algorithm 

Quantum decoherence was incorporated into the simulation through the Consistent His- 
tories theory discussed above. The computational algorithm given below follows directly 
from Eq. |Tl] and the coherence time intervals were chosen from a Poisson distribution with 
characteristic time scale r/5 = 3.1 fs which was derived from our previous estimates of the 
coherence time scale for an electron in H2O when the excited state is nearly solvated. As 
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discussed in a previous paper, this time scale was estimated by computing the average de- 
cay time of the overlap of frozen Gaussian vibrational wavefunctions evolving on different 
adiabatic potential energy surfaces by sampling a large number of nearly equivalent excited 
state configurations. This is in marked contrast to previous simulations on this system in 
which coherences between electronic states were not maintained beyond 1 fs. The consistent 
histories algorithm proceeds as follows: 

1. Determine new coherence interval from Poisson distribution of possible intervals with 
characteristic time scale T£,. 

2. Propagate the quantum wavefunction over this time scale while self consistently evolv- 
ing the classical degrees of freedom. At the end of each dynamical time step h < t^, 
determine the switching path followed by the classical variables using the stationary 
phase algorithm developed by Webster, et a/.0il modified such that coherence in the 
quantum wavefunction is maintained over the entire coherence interval. We also note 
that the stationary phase switching path is determined in a piece-wise continuous 
fashion by selecting intermediate quantum states every dynamical time step. This is 
to avoid the computational overhead of computing variationally the stationary phase 
trajectory over relatively long time intervals. So long as To is no longer than a few 
dynamical time steps, this approximation should not be too extreme; however, certain 
pathological cases can be invented in which this approximation does break down. 

3. If either a switch occurs in the time interval or we reach the end of the interval, the 
quantum wavefunction is collapsed using the projection operators discussed above. 

4. Repeat. 

B. Results 



In Fig. we plot the switching times from the excited state to the ground state for 
a total of 23 simulation runs. The starting configurations were generated by performing 
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a 35 ps simulation in which the electron was prepared in the ground state and 16 initial 
configurations were chosen whenever the energy diff'erence between the ground state and 
one of the p-like excited states become resonant with the excitation laser (1.7 eV). llUll 
These starting configurations were identical to configurations used previously by Schwartz 
and Rossky in their work on this system. llHil Following the initial excitation, the system 
was allowed to evolve. During this time, the energy gap between the p-states and the 
s-state narrowed to ~ 0.5eV as the excited state was solvated by the surrounding water 
molecules. The simulation continued until a switch from the excited state to the ground 
state was recorded. Immediately after the switch, the energy gap between the occupied 
ground state and the first excited state widened dramatically as the solvent responded to 
the new electronic state. The average and median switching times from these simulations 
are compiled in Table 1. 

In order to test the sensitivity of the switching times to the choice of sequence of coherence 
intervals, five configurations were "recycled" by using different random number sequences 
for the coherence time intervals to generate different switching paths starting from the same 
initial configuration. For each such pairs of paths, the classical dynamics were identical 
up until the earlier switching time when one of the paths switched from the excited state 
to the ground state. However, since the coherences between the different quantum states 
were killed off at randomly different times for each path, each path sampled a different 
probability distribution function for making a switch at each MD time step. Hence any 
correlations between switching times resulting using the same initial configuration refiect 
correlations in the distribution of coherence time intervals. Given that there is typically a 
100 fs time difference between pairs of data (in one case a 530 fs difference) which is roughly 
1/3 of the average survival time scale of the excited state, we find very little correlation 
between the switching times originating from the same initial configuration. 

The excited state survival probability refiects not only the probability of making a tran- 
sition from the excited state to the ground state, it also refiects the solvation dynamics of 
the non-equilibrium excited state following its preparation. Initially, the energy gap between 
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the excited and ground state is very wide and a non-adiabatic transition between the states 
is highly unhkely. However at later times, after the excited state has begun to be solvated 
by the surrounding water molecules, the energy gap closes considerably and nonadiabatic 
transitions become more likely. The solvation dynamics of this system have been extensively 
studied previously by Schwartz and Rossky.iiil^ Their results for the solvation response, 
defined as the normalized autocorrelation of the energy gap between the ground and excited 
states, indicates a rapid 24 fs initial Gaussian component which makes up roughly 40% of 
the response followed by a longer time 240 fs exponential decay which makes up most of the 
remainder. 

As the excited state is solvated by the surrounding water molecules, the energy gap 
between the excited state and the ground state narrows to its equilibrium value. Further- 
more, from first order perturbation theory, we expect that the electronic transition rate to 
be inversely proportional to the magnitude of the energy difference between initial and final 
states. Thus, a simple model for the excited state survival probability can be written 



Pit) (16) 



dt iu{t)' 

subject to the initial condition P(0) = 1. Here, kgq is the non-adiabatic transition rate 
for the solvated excited state and uj{t) is the average energy gap normalized to the average 
energy gap of the solvated excited state. U{t) is also related to the solvation response S{t) 
by 

^^^^ ^ (^o)g(t) + Hq)(l-g(t)) _ ^^^^ 

The solvation response can be modeled as short time Gaussian decay followed by a longer 
time exponential relaxation,&B0 

Sit) = Ae-i/2(tAG)^ + (1 _ A)e-'/'^. (18) 

In Fig. ^, we plot uj{t) using a parameteric fit to the solvation response function obtained 
from our simulations: A = 0.38, to = 0.38 fs, tE = 240fs. These parameters agree exactly 
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with those obtained by Schwartz and RosskycJiHa for this system. This comes as no surprise 
since our inclusion of quantum coherence effects does not effect the computation of the 
interaction forces between the excited electron and the water molecules. At short times, 



At long times Lj(t) — 1 as the energy gap relaxes to its equilibrium value. Thus, k^g is the 
exponential decay constant of the excited state population once the excited state is fully 
solvated. 

The solution of Eq. [l^ is the one parameter family of curves given by 



Using a non-linear fitting procedure, we fit our data to this model to obtain an estimate 
of the equilibrium nonadiabatic lifetime of 241 fs with a = 1-707. A plot of our data 



simple model as nicely as the data given by Schwartz and Rossky.cJiHScZl This suggests that 
the longer coherence times used in this study imparts a non-trivial memory dependency into 
the survival probability which is inadequately captured in Eq. 

As expected, the lifetimes reported here (Table 1) are consistently shorter than the life- 
times reported by Schwartz and Rosskyi^hBS in which a constant 1 fs coherence time scale 
throughout their simulations, thus emphasizing the profound sensitivity of these simulations 
to the coherence time scale. Furthermore, our results are consistent with the estimated 
lifetimes reported in SBPR@ where the effective equilibrium lifetime was estimated as a 
function of coherence time by reanalyzing the transition amplitudes computed in the origi- 
nal Schwartz and Rossky simulations with the coherence time scale estimated from a frozen 
Gaussian approximation to the nuclear wavefunction. 



when U{t) > 1, the nonadiabatic transition rate will be small since the energy gap is large. 




(19) 




IV. DISCUSSION 



In this paper we have briefly described the results of our work towards a molecular level 
description of quantum relaxation phenomena. Here we have focused exclusively upon the 
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role that transient quantum mechanical coherences between the solvent and the solute play in 
the electronic relaxation of an excited solute species. In mixed quantum-classical computer 
simulations, fundamental assumptions about the decay of these coherences produce direct 
manifestations on the computed quantum mechanical transition rates and must be included 
consistently in order to make realistic predictions and comparisons.i'0 

The computational algorithm, based upon the so-called consistent histories interpreta- 
tion of quantum mechanics, provides both the molecular level underpinnings of quantum 
decoherence and the computational means for properly including decoherence effects in 
non-adiabatic quantum-molecular dynamics simulations. According to the rules of ordinary 
quantum mechanics, a quantum system will evolve into a coherent superposition of alter- 
native states. In our decoherence theory, this coherence is dissipated due to the differences 
in the forces exerted on the bath by alternative states involved in the superposition. Thus, 
on short scales, a quantum system in a bath will obey the rules of ordinary quantum me- 
chanics and evolve into a coherent superposition of states, whereas on longer time scales, 
the coherences between states are diminished and the quantum system must be described as 
statistical (i.e. incoherent) mixture of states. As the quantum system interacts continuously 
with the bath, coherences between states are continuously created by non-adiabatic coupling 
and damped by the divergence in the bath dynamics induced by the system-bath coupling. 

This subtle interplay between coupling and decoherence and the subsequent dependency 
of transition rates on the decoherence time scale has profound implications for a variety of 
condensed phase chemical dynamics including: internal conversion and internal vibrational 
energy redistribution (in which the bath is comprised of all the modes of the molecules except 
for the one mode of interest), electronic energy transfer between molecules or different parts 
of the same molecule, and charge transfer reactions including proton and electron transfer. 
In these latter examples, both the condensed phase environment and the internal motions 
of the molecules act as a bath which couples the quantum states together. The decay 
of quantum coherence, which depends upon the frequencies and populations of the bath 
modes coupled to the quantum system will determine the extent to which the non-adiabatic 
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coupling can act to allow the chemical reaction to proceed. Changes in the spectral density 
due to isotopic changes in the bath can have a substantial impact on non-adiabatic chemical 
dynamics.i Furthermore, the decay of quantum coherence can determine the adiabaticity 
for a chemical reaction. 

Perhaps the two major lacunae in our present theory is the explicit dependency upon 
an a priori estimate of the coherence time scale and the fact that this time scale remains 
fixed throughout the calculation. As mentioned above, we estimated this time scale by 
computing the average decay time of the overlap of a product of Gaussian coherent states 
evolving on different adiabatic potential energy surfaces. In this estimate, the individual 
widths of the coherent state wavefunctions centered about the initial phase space points of 
the classical nuclei are set to be proportional to the thermal DeBroglie wavelength of each 
nuclei. Although physically realistic for a variety of situations, this does leave the coherence 
length scale (i.e. the widths) as an adjustable parameter. While the effects of changing the 
coherence length scale over a broad range have not been systematically studied, results from 
our previous work demonstrate that the non-adiabatic transition is quite sensitive to changes 
in the coherence time scale and hence will be sensitive to changes in the coherence length 
scales. Furthermore, as the bath explores various regions of the quantum potential energy 
surface, the force differences which give rise to the decay of the quantum coherences^ should 
vary from one configuration to the next. Current progress is underway towards obtaining 
both the coherence length and time scale during the course of the non-adiabatic simulation 
ab initio by examining the quantum mechanical fluctuations of the bath particles about their 
stationary phase paths. 
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Tables 



TABLE I. Excited state lifetimes for the aqueous electron. (SR= Schwartz and Rossky, J. Chem. 
Phys. 101 6902 (1994). SBPR = Schwartz, Bittner, Prezhdo, and Rossky, J. Chem. Phys. 104 
4942 (1996)). The coherence times used in each study are listed in parentheses. 





Proscuit (3.1 fs ) 


SR (1 fs) 


SBPR (2.8-3.1 rs)('') 


Median 


338 fs 


630 




Average 


384 


730 




E(nulil)rium 


231 


150 


310-270 



(a.) Only Equilibrium lifetimes considered. 
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Figure Captions 

FIG. 1. Energy eigenvalues and quantum populations along different switching paths for the 
electron in water. Shown in the upper graph are the energy eigenvales of the occupied electronic 
state for two possible switching paths with the occupations shown in the lower plot. Following 
promotion to the excited state, the electron makes a series of stochastic hopping attempts between 
the initial and final states with two possible outcomes shown. 

FIG. 2. Excited state cumulative survival probability as a function of time following initial 
excitation. Superimposed points (A) are the switching times from the simulations. Curve B is the 
fit of the data to the energy gap relaxation model given by Eq. 22 and curve C is a fit of the data 
to a Gaussian. See text for details. 

FIG. 3. Plot of the normalized energy gap, oo, between the ground and occupied excited state 
following initial excitation. 
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